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ABSTRACT 

Context. One of the greatest difficulties in astrophysics is the determination of the fundamental stellar parameters, one of which is 
the initial mass fraction of helium (Y). However, given that Y can be measured spectroscopically in only a small percentage of stars, 
a linear relationship is assumed between Y and the mass fraction of metals (Z) from a canonical perspective of the chemical evolution 
of the galaxies. This Y-Z relation is generally represented as Y = Yp + AY/AZ x Z, with the value of the helium-to-metal enrichment 
ratio (AY/AZ) assumed as a constant. However, there is no fundamental reason for eveiy star to have a Y value on a linear scale with 
Z. Indeed, different AY/AZ values may be expected in different populations which have undergone different chemical enrichment 
histories. 

Aims. In this paper a new method for determining the fundamental stellar parameters of nearby stars is presented that uses at the same 
time Mboi, T^ff, and log g. One of these parameters is Y, which is used to determine the validity of the Y-Z relation. 
Methods. A new set of evolutionary tracks is created using the PGPUC stellar evolution code, which includes 7 masses (0.5 < 
M/Mq < 1.1), 7 helium abundances (0.230 <Y < 0.370), and 12 metallicities (1.6 x lO""* < Z < 6.0 x 10"^) for solar-scaled chemical 
compositions ([ff/Fe] = 0.0). The suggested method is tested using two different spectroscopic databases of nearby main sequence 
stars with precise parallaxes, and spectroscopic measurements of [Fe/H], Teff and g. 

Results. The proposed method is compared to other techniques used to determine the fundamental stellar parameters, where one 
assumes an age of 5 Gyr for all nearby stars. This comparison demonstrates that the hypothesis regarding constant age leads to an 
underestimation of the Y value, especially for low metallicities. In addition, the suggested method is limited to masses above 0.60 
Mq and requires high-precision measurements of spectroscopic surface gravities in order to obtain reliable results. Finally, estimating 
masses and ages assuming a Y-Z relation rather than a free Y value may induce average errors of approximately 0.02 Mq and 2 Gyr, 
respectively. 
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1. Introduction In its simplest (and most widely used) form (e.g. . 
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Peimbert & Torres- 

The relationship between the mass fraction of helium Y and IWilson & Rood||l994 andreferences therein), the Y-Z relation 

the mass fraction of metals Z, hereafter referred to simply reads as 7 = Yj, + AY/AZ x Z, wher e Y„ is th e primordial he- 

. ^ as the Y-Z relation, is of paramount importance in studies Hum abunda nce (e.g., Yp * 0.240; Izotov et al.||2007; Steigma^ 

H of the formation and evolution of most of the visible com- |2007i |2012t and references therein), and AY/AZ is the helium- 

•^.ponents of the Universe, including planets, stars, star clus- to-metal enrichment ratio. This relationship can also be written 

ters, and galaxies. That stellar evolution depends strongly as Y = Yq + AY/AZ x (Z - Zq), where Yq and Zq are the solar 

on Y has been well established since the 1960s, when stel- mass fraction of helium and metals, respectively, 

lar structure and evolution calculations revealed that main se- The measurement of helium abundance in stars at different 

quence stars with high helium abundances are brighter, hotter, metalUcities is a direct method for calibrating the Y-Z relation, 

and evolve more rapidly than their low-Y counterparts (e.g.. However, given that strong helium lines are only present in stars 

Demarque 1967; Iben & Fatilk iier 196 8; Si moda & Iben 1968, with effective temperatures (reft) higher than 8 000 K, which are 

1970l; jPem argue et al. 1 19711: iHartwick & VandenBerd 119731; only reached by main-sequence stars with masses greater than 



Sweigart & Grossiil978i) . One of the most important, though of- a; 1 .5 Mq, the chance of finding sta rs with these temperatures at 



ten unstated, appUcations of the Y-Z relation is in relation to low metallicities is low dSuda et al ll200i) . since low-metallicity 

theoretical models, in which a free parameter (the helium abun- stars evolve faster (e.g., Simoda & Ibenl ll968l : lAizenman et al.l 

dance) is avoided and it is assumed that divergences between 1969; Iben 1974; Simoda & Iben 1970). Exceptions are the blue 

theory and observations are only associated to differences in age, horizontal branch (BHB) stars present in some globular clusters, 

mass, and/or other free parameters, such as the metallicity. which are old, metal-poor stars in the central helium burning 
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phase. However, BHB stars with T^ff higher than 1 1 500 K show 
higher metal and lower helium abundances (effects that are inten- 
sified for higher Teff, Grundahl et al. 1999; Behr 2003), probably 
as a consequence of diffusion/levitation of elements. This leaves 
a narrow band of T^s in which the initial helium abundance of 
these stars can be measured, although this remains complex due 
to the weakness of the lines and the high S/N required. This has 
been demonstrated by iVillanova et al. (2009bi) in stars of NGC 
6752, employing the He i line at 5876 A. Using the same He line, 
IVillanova et 51(1200931) determined a high overabundance of he- 
lium with respect to the Sun in the hottest main sequence stars 
of the open cluster NGC 6475 ([Fe/H] = -1-0.03). Alternatively 
based on the chromospheric Hei line at 10830 A. iDupree et al.l 
(1201 Ih established the existence of variations in the helium abun- 
dances of red giant branch (RGB) stars of oj Cen with similar 
properties (Teff, brightness, and metallicity). Almost simultane- 
ously and with the same line at 10830 A. iPasquini et alJ (1201 ih 
obtained similar results comparing two RGB stars in the glob- 
ular cluster NGC 2808. Despite these significant efforts to de- 
termine helium abundances spectroscopically, present-day high- 
precision observations are not sufficiently accurate to study this 
element in detail for a large amount of stars. 

Since direct He measurements are not yet available for most 
stars, indirect methods have been used to calibrate the Y-Z re- 
lation semi-empirically. One such calibration is performed us- 
ing a standard solar model and by assuming that the Sun rep- 
resents all stars with solar metallicity. With this in mind, by 
applying the value of Yp for Z - and the solar calibration 
of Yq and Zq, one can obtain the AF/AZ value. In o ur case, 
Yq = 0.262 and Zq = 0.0167 dValcarce et al.ll2012l) . giving 
the result AF/AZ = 1.31. However, Yq and Zq depend on the 
input physics of the theo retical model (Pietrinfemi et al. 2004,; 
VandenBera et al. 2006; IWeiss & Schlattri2008l; iBertelh et all 
2008i; iDotteret al., ,20081) and the a s sumed solar chemical dis- 
tribution (|Grevesse &Noelsl 119931: iGrevesse & Sauvall Il998l 



lAsplund et alJl2005Ll2009HCaff^au et alJl201 ih 



Another sem i-empirical determination of the Y-Z relation 
is carried out bv ICasagrande et aL (|2007[ hereafter C 07) using 
the sample of K dwarf stars in Casagrande et al.l (l2006i hereafter 
C06) with known absolute bolometric magnitudes (Mboi). This 
method involves obtaining helium abundances using isochrones 
in the theoretical Teff - Mboi plane. Assuming that all the stars 
in their sample are 5 Gyr old, the authors determine AF/AZ - 
2. 1 ± 0.9 for stars with solar and above-solar metallicities. How- 
ever, the helium abundances estimated in metal-poor stars by 
these authors are too low to be considered real. U sing a simi- 
lar method with K dwarf stars iGennaro et al" l (l20T0l) determine a 
AF/AZ = 5.3 + 1 .4, concluding that the assumption of a constant 
age for all stars in the sample leads to underestimation of AF/AZ. 
This is because stars evolve faster for high helium abundances 
and/or low metallicities. Other studies usin g similar methods to 
determine A F/AZ include iFaulknej (1 19671) with AF/AZ = 3.5, 



|Penin et al] (Il977l) with AF/AZ = 5 . iFernandes etaP (1 19961) 



wit h AF/AZ > 2.|Pagel & Portinaril(ll998l) with AF/AZ =3+2 



and|Jimenez.etalJ (|2003|) with AF/AZ = 2.1 + 0.9 



An add i tional method to determine AF/AZ is that presented 
bv lRenzinil (1 1 9941) . who determines 2 < AF/AZ < 3 using the ra- 
tio of horizontal branch clump stars to red giant stars in the Milky 
Way's bulge, compared to the theoretical predictio n for different 
F. Using a similar approach, ISalaris et al.l (l2004l) determine an 
initial helium abundance of F = 0.250 + 0.006 for 57 globular 
clusters in a wide range of metallicities (-2.2 < [Fe/H] < -0.3), 
that is, AF/AZ * for globular clusters. 



Of course there are other methods for determining AF/AZ. A 
comprehensive review of these different methods can be found in 
Sect. 9 of Gennaro et al. (2010). The different available methods 
can be summarized as follows; i) predictions of galactic chemi- 
cal evolution models (1.5 < AF/AZ < 2.4), ii) analysis of the 
chemical composition of planetary nebulae (2.0 < AF/AZ < 
6.3), and iii) measurement of helium recombination lines of 
Galactic and extragalactic Hii regions (1.1 < AF/AZ < 5.0). 
However, taking into account only the last ten years, there is a 
general consensus favoring a AF/AZ value around a 2.0. 

In this paper a method for determining AF/AZ and the main 
fundamental stellar parameters of nearby field stars is presented. 
Section |2] details the theoretical evolutionary tracks created for 
this purpose, which are used to obtain AF/AZ (Sect. O. Sec- 
tion|4]briefly demonstrates an application of this method for two 
databases of nearby stars, with conclusions presented in Sect. |5] 



2. Stellar models 

The stellar models used in this paper are calculated using the 
PGPUC stellar evolution code dValcarce et al. 11201 2i hereafter 
PGPUC SEC) that was recently updated with the following 
input physics: radiative opaci ties for high (llglesias & RogersI 
1 19961) and lo w temperatures ([Ferguson et al. 120051) : conduc- 
tive opacities ('Cass isi et alj 20071) ; thermonuclear reaction rates 



(Anguloetal. 19991; iKunz et alj 12002; Formicola et al. 2004; 



Imbriani et al. 2005 ); equ ation of state (Irwin 20(j7l) ; mass loss 
dSchroder & Cuntz, |2005|) ; and boundary conditions (ICatelanI 



The PGPUC SEC was also calibrated using the Sun for a 
stellar model with a sola r mass and a solar chemica l composition 
distribution according to lGrevesse & Sauvall (1 19981) to reproduce 
the prese nt solar luminosity, rad ius, and {Z/X)q - 0.0231 ±0.005 
ratio (G revesse & Sauvailll998l) . The results indicate a mixing 
length parameter »/ - 1.896, an initial solar helium abundance 
Fq = 0.262, and a global sol ar metallicity Zp = .0167. For a 
more detailed description see IValcarce et al.l (1201 2lP I. 

A set of stellar models was calculated using a solar-scaled 
distribution of alpha-elements ([a/Fe] - 0.0). These stellar 
models consider the evolution of stars from the zero-age main 
sequence (ZAMS) to the tip of the red giant branch. The 
grid of initial parameters includes i) 7 masses, between 0.5 
and 1.1 Mq at intervals of AM - 0.1 Mq; ii) 7 helium abun- 
dances, for F = 0.230 and from F = 0.245 to F = 0.370, 
at intervals of AF = 0.025; and iii) 12 global metallicities 
(Z = 0.00016, 0.00028, 0.00051, 0.00093, 0.00160, 0.00280, 
0.00503, 0.0089, 0.01570, 0.01666, 0.03000, and 0.06000). For 
this range of Z, iron abundances with respect to the Sun change 
from -2.05 < [Fe/H] < 0.56 when assuming F = 0.245, to 
-1.96 < [Fe/H] < 0.66 for F = 0.370. This variation is im- 
portant for a reliable determination of the stellar parameters (see 
next section). 

The present study also used the set of stellar models with 
an enhanc ed distribution of alph a-elements [a/Fe] = -1-0.3 cal- 
culated in IValcarce et al.l d2012l) . The grid of properties in- 
cludes the same initial masses and helium abundances listed 
for [a/Fe] - 0.0, and a set of 9 global metallicities covering 
1.6x10-'* <Z< 1.57x10-2. 



' All the stellar models calculated for this paper, or those interpolated 
from these models (for any mass, ¥, Z, or [a/Fe]), can be downloaded 
from .http : / /www2 . astro . puc ■ cl/pgpuc i 
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The zero-point for Mboi wa s determined using the respective 
solar value Mboi.0 = 4.77 mag jBahcall et alj|1995|pl 

3. Semi-empirical determination of fundamental 
stellar parameters 

The fundamental stellar parameters (FSPs) are all the parame- 
ters of a star defining its present-day observational properties. 
The main FSPs include the initial mass (M), the initial chemi- 
cal composition (usually Z, 7, and [a/Fe]), and Age. However, 
other parameters can also affect the star's present-day observa- 
tional properties, including the rotational velocity, convective 
overshooting, magnetic field, mixing length parameter, and mass 
loss rate (after the main-sequence phase), among others. With 
this in mind, establishing the FSPs using theoretical models con- 
strained by observational measurements is not a straightforward 
task. This is because the number of FSPs is generally greater 
than the amount of observational properties that can be mea- 
sured, resulting in a degenerate mathematical problem. Thus, 
the present study aims to select the most important FSPs that 
match the number of observational constraints available. 

In online databases, common observational properties of 
nearby stars with available parallaxes are i) Mboi, or the absolute 
magnitude in a given filter; ii) Teff ; iii) the spectroscopic surface 
gravity (gfspec); iv) the iron abundance with respect to the Sun 
([Fe/H]); and in some cases v) the distribution of alpha-elements 
with respect to the Sun ([ff/Fe]). Using these observational prop- 
erties, this paper aims to estimate M, Age, Y, and Z for each star, 
and then obtain ISYIhZ. When [Fe/H] and [a/Fe] are known, it 
is possible to restrict the theoretical relationship between Z and 
ras 



[Fe/H], ^[M/H]+ log 



fFe(O) ^ mzia) 



(1) 



where [M/H] = log ( i-F-z )~^°§(f ^ve{a) is the number frac- 
tion of iron with respect to all the elements heavier than helium; 
mzia) is the average atomic mass of heavy elements weighted 
by the number of atoms, and {Z/X)q is the solar ratio of metals 
with respect to hydrogen. It is important to note that fFe(a') and 
mzia) depend on [a/Fe], and fFe(O) and mz{0) correspond to the 
case [a/Fe] - 0.0. The values of the parameters that depend on 
[a/Fe] are shown in Table [T] 

Based on the adopted reference mix, one can also obtain an 
expression relating the overall quantity of metals (as commonly 
represented by Z or [M/H]) to that contained in the form of a- 
capture elements and iron only (as commonly expressed in terms 
of [a/Fe] and [Fe/H], respectively). In other words, we seek to 
derive the a and b coefficients for which the following expression 
is valid, for the adopted chemical composition 



Zoj lZ„ J IZPe. 



(2) 



^ A more recent Mboi.o value is 4.7554 ± 0.0004 mag (see Eric Ma- 
majek's web page http://www.pas.rochester.edu/~emamajek/sun.txtl, 
which does not affect the main results, thanks to the small difference 
with the value considered in this work. 



-[Fe/H] = -0.17 [a/Fe] = 0.0 

_ Y Z Mass Age 

[M„] [Gyr] 
-0.230 0.0118 0,844 0,0 

0.245 0,0116 0,819 0,1 
-0,270 0,0112 717 15,1 

0,295 0,0108 0-636 29,3 
-0.320 0.0105 0.573 42.5 

0.345 0.0101 0.517 55.5 
-0.370 0.0097 0.500 47.1 




ZAMS 
Y-0.230 
Z.0.0118 





3.76 3.74 3.72 3.7 
log T.„ 



3. 68 



4.3 4,4 4,5 4.6 
log E 



Fig. 1. Method for determining the FSPs, as applied in this case to a 
star with (log T^g, Mf,a])siar = (3.7, 6.0). Left panel: First, for each star 
with a given [Fe/H], the masses and ages of evolutionary tracks (contin- 
uous lines) are determined for a group of helium abundances (maintain- 
ing [Fe/H] constant) at the desired point. Upper-right panel: Zoom-in 
of the left panel. Next, the helium abundances with reliable matches 
are selected. In this case, two Y values are disregarded i) Y = 0.230, 
because the star is fainter than the ZAMS for this helium abundance 
(dotted line), and ii) Y = 0.370, since a mass lower than 0.5 Mq (the 
lower mass limit of the set of tracks) is required to pass through the 
point (log Teff, MboOiror- Lower-right panel: Since Y displays an al- 
most linear dependence with log g at the given (log Teff . ^boi).«nr point, 
Y can be determined using the spectroscopic g value. The continuous 
and dotted vertical lines represent the spectroscopic log g value of the 
test star and its error, respectively. Dots represent the theoretical log g 
value at the point (logT^a, Mboi),^,. for each Y value, and their errors 
are obtained assuming errors in Tcs and Mboi of CTr^g = ±50 K and 
iTMf,^! = ± 0.05 mag, respectively. 

very si mila r to the one that was found previously bv lSalaris et al.l 
(il993h and iYi et alJ (1200 Ih . namely 



[M/H] = [Fe/H] + log (0.6355 x /„ + 0.3646), 



(3) 



with/„ = lOt'^/''"]. 

The first step of our method is to determine the relationship 
between Z and Y for each star with a given [Fe/H] and, indirectly, 
a given [a/Fe]. From Eq. |3]it is straightforward to obtain the 
following relationship 



with a+b - 1, and where Z, and Zpe represent the mass fractions 
of a-elements and iron, respectively. Assuming that O, Ne, Mg, with 
Si, S, Ca, and Ti are all a-capture elements, we find a - 0.6355 
and b - 0.3646. This leads naturally to an expression that is 



1 +c 



c = (0.6355 X /„ -I- 0.3646) x (-) x 10^^"^"'. 



(4) 



(5) 



The second step is to determine the other FSPs depending 
on Y: Z{Y), M{Y), and hgtiY). To that end, a set of interpo- 
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Table 1. Parameters used to estimate [Fe/H]. 



[a/Fe] 


-0.2 


+0.0 


+0.2 


+0.4 


+0.6 


+0.8 




0.031012 


0.023495 


0.016988 


0.011787 


0.007936 


0.005229 


mz{a)/mz{0) 


0.989131 


1.000000 


1.009287 


1.016946 


1.022616 


1.026603 


o L tFc(O) mz(a) J 


0.115808 


0.000000 


-0.136819 


-0.292274 


-0.461661 


-0.641154 



lated evolutionary tracks is used in the theoretical Hertzsprung- 
Russell diagram. As is shown in Fig. [T] the stellar evolution- 
ary tracks with different Y values pass through the test star at 
log Tetf = 3.7 and Mboi = 6.0. However, as can be seen, there are 
evolutionary tracks that do not pass exactly through this point for 
some Y values, given that i) the star is fainter than the ZAMS for 

Y = 0.230 (dotted line in the upperright panel ofFig.[T]i, indicat- 
ing there is no mass value to solve the problem for this specific Y 
value; or ii) the star is cooler (hotter) than the whole evolution- 
ary track with the minimum (maximum) mass of our database, 
as is the case for Y - 0.370 in that figure. Note that an evolu- 
tionary track with a mass lower than 0.5 Mq for Y = 0.370 can 
pass through the point (log T^ff, Mbni)siar, but because stars with 
masses lower than our limit are more difficult to detect and have 
large associated observational errors, these evolutionary tracks 
are not calculated. Note also that these stars are much more 
difficult to model thanks to the molecules formed in their cool 
atmospheres. In both cases, the FSPs for these specific Y values 
are ignored becuase of physical inconsistency. 

Finally, since each evolutionary track with a different Y value 
has a different mass for the same log T^s-MboX point, the helium 
abundance can be determined by comparing gfspec with the theo- 
retical surface gravity {g). As is known, in the Mboi-T^eff diagram 
the differences in g are only due to the different stellar masses, 
since the stellar radius is determined by the Stefan-Boltzmann 
law. However, when a specific filter is used instead of Mboi 
(for example My), the relationship is not entirely the same since 
bolometric corrections do not have a linear dependency with T^s 
and g. 

The bottom-right panel of Fig. [T] shows the relationship be- 
tween Y and g with the same Mboi and Teff values, where each 
point is the value of log g at the test point of the track with a 
given Y value. Note the almost linear relationship between Y 
and log g. Theoretical errors for log g are determined assuming 
typical errors in T^s and Mboi of ±50 K and ±0.05 mag, respec- 
tivel y, based on present-da y measurements of these parameters 
(e.g.. lBaumann et alj2010h . Using the spectroscopic value of the 
test star log g - 4.53 ± 0.06 dex, values represented by vertical 
continuous and do tted l i nes, a nd the Hermite interpolation algo- 
rithm presented by [Hill (UMll); the following FSPs are obtained: 

Y = 0.287!00^i. z = 0.0109!0°ojj5. ^ ^ 0.661!»igMo; and 

Age = 24.7!;^;^ Gyr 

Since the separation of evolutionary tracks in the Mboi-rgif 
diagram at the red giant branch is far smaller than the typical 
errors, this method must be restricted to main-sequence and sub- 
giant branch stars, although the latter may have a larger uncer- 
tainty in the derived age. As such, stars of the databases men- 
tioned in the following section are restricted by log g > 3.8 dex. 

3. 1 . Comparison of methods 

In this section the method described above is tested and com- 
pared to other similar approaches. The results are shown in Fig. 
12] where the solar symbol and the dotted line represent the Sun's 



properties and the Y-Z relation Y - 0.240 + 2.0 x Z, respec- 
tively. Extrapolated results are used in this section only; the risks 
involved in this procedure are discussed in Sect. 13.21 

Stars listed in the C06 catalog are used, consisting of nearby 
low-mass main sequence stars with good [Fe/H] and T^ff mea- 
surements. However, since the method suggested in this paper 
uses the ^spec value, the C06 table had to be completed with 
the spectroscopic values obtained by the author of the refer- 
ence cited for each star The range of iron abundances, effec- 
tive temperatures, and surface gravities covered by this database 
is -2.0 < [Fe/H] < +0.4, 4400 < reff[K] < 6400, and 
4.1 < logg < 5.0, respectively, with errors of approximately 
C[Fe/H] ~ ±0.15 dex, cTj-.ij. < ±100 K, and o-\ogg < ±0.20 dex, as 
well as some outliers. For stars without an estimated log g error, 
a criogg = 0.05 dex is assumed. The apparent bolometric magni- 
tude mboi calculated by C06 and the parallaxes for each star are 
also used to determine the absolute bolometric magnitude Mboi. 

The first method is similar to that described by C07, used to 
obtain the FSPs of the C06 stars. As previously mentioned, these 
authors assume that all stars are 5 Gyr old, and determine the 
helium abundance for each star using isochrones in the theoret- 
ical Mboi-7"eif plane. Given that only a limited set of isochrones 
with different helium abundances is available in the literature, 
C07 use extrapolated results to obtain Y values between 0.10 
and 0.30. In this first approach (upper panels of Fig. |2), rather 
than obtaining the FSPs using g^pec (as described in the previous 
section), here the Age is used as the known parameter to inter- 
polate (or extrapolate) the other FSPs for a fixed value of 5 Gyr. 
The results indicate that almost all stars with Z > 0.01 follow 
the adopted Y-Z relation. However, for Z < 0.01 the Y value 
is significantly dispersed, with a trend towards lower Y values 
for lower metallicities. This trend is also found by C07, but they 
attribute it to the current limits in stellar models. 

There is an alternative explanation for the steeper trend be- 
tween Y and Z at low metallicities. Stars that are older than 
5 Gyr will be brighter than theoretical models with that age 
and, consequently, mistakenly appear He-depleted. This is 
because when Y decreases the main s equenc e locus becomes 
cooler and seems brig h ter (R emarque] Il96"7l ITb en & Faulkne^ 
'1968*; 'Simoda & Iben 1968, 1970; Demaraue et all 1197 iT 
iHartwick & VandenBerg. ,1973: Sweigart & Gross, ,19781) . Al- 
though this effect is not so important for stars with masses lower 
than 0.60 MeB when the mass is greater than 0.60 Mq the evolu- 
tion is faster and the underestimation of helium is higher Sim- 
ilarly, stars younger than 5 Gyr will appear overabundant in he- 
lium. These effects are amplified when met allicity decreases or 
helium increases (see IValcarce et al.l I20T2I for a recent pub- 
lished explanation). As shown in the upper-right panel of Fig. |2] 

^ Stars with masses lower than 0.60 Mq show almost no M],„\ and T^g 
variations, among other properties, over a period equivalent to the age 
of the Universe. Using PGPUC SEC calculations for a 0.60 Mq star 
withZ = 0.01, Y = 0.245, and [a/Fe] = 0.0, the variations in Mboi, T^ff, 
and g after 13.5 Gyr are A(Mboi) a 0.16 mag, A(logre(f) ~ 0.01 dex, 
and A(log g) s; 0.03 dex, respectively. 
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almost all stars have masses greater than 0.60 M©, and accord- 
ingly the ages of nearby stars cannot be assumed similar to the 
Sun's age for comparisons of this kind. 

In the middle panels of Fig. |2]the FSPs of stars using the 
method proposed in the previous section are shown. In this case 
it is assumed that the a-element distribution of each star was not 
determined, since only a small fraction of spectroscopic studies 
have provided estimates of [ff/Fe] for each star. Thus, the FSPs 
are obtained using evolutionary tracks with a solar-scaled distri- 
bution of elements ([o'/Fe] = 0.0). In this panel, stars older and 
younger than 13.5 Gyr are represented by circles and triangles, 
respectively. Results obtained by interpolation and extrapolation 
are depicted by filled and open symbols, respectively. The most 
notable difference from the previous approach (all stars are as- 
sumed to be 5 Gyr old) is the reduced number of stars with very 
low helium abundances. Moreover, there is a high percentage 
(68%) of stars older than 13.5 Gyr, which is far more evident 
for lower masses. There are two reasons for this: i) the cur- 
rent problem in predicting reaUstic radii for low-mass stars (e.g. , 
iTorres et al.ll2010t iFeiden & Chabovetjl2012l iBasu et al.ll2012l) 
where errors in the radius determination of only 5% produce dif- 
ferences of some tens of Gyrs, and ii) any small error in the 
estimations of Mboi, Teft, and/or g can induce a large error in the 
age, due to the small variations expected for low-mass stars dur- 
ing a period equivalent to the age of the Universe. These prob- 
lems induce a lower mass limit of ~ 0.60 for the reliability 
of estimations of Y and Age with the suggested method because 
of the high sensitivity of these parameters to small variations of 
^boh Teff, and/or g. Alternatively, masses are less sensitive to 
variations of these observational properties, thus making mass 
estimations more reliable. 

The bottom panels of Fig. [2] show the results of the most 
robust method, which requires the spectroscopic determination 
of [a/Fe] for each star. This method applies the same procedure 
described in Sect. |3] with a set of evolutionary tracks according 
to the [o'/Fe] value of the respective star. This set of evolution- 
ary tracks for each [a/Fe] is interpolat ed in accordance with the 
procedure described in the appendix of lValcarce et al.l(l2012h . In 
order to have the same number of stars as in the middle pan- 
els, stars with [cc/Fe] < 0.0, and [c/Fe] > -hO.3 are assumed to 
have [ff/Fe] values of 0.0, and -t-0.3, respectively. As demon- 
strated, when the [a/Fe] value is taken into account, dispersion 
on the initial helium abundance decreases even further than in 
the previous case, particularly for low metallicities. Although 
this method is more robust than the previous one, results still in- 
dicate ages greater than the age of the Universe for stars with 
masses M < 0.70 Mq. Moreover, the statistical difference be- 
tween the two methods (considering or not [c/Fe]) is not signif- 
icant, at least for this set of stars. 

3.2. Why extrapolated results should not be used 

Although in some mathematical problems results can be ob- 
tained using extrapolation algorithms, these results may not al- 
ways be realistic. For example, when the FSPs of stars are deter- 
mined using their log g^pec value (as shown in Fig. [Ij, there are 
cases when gfspec values are greater than g of the minimum avail- 
able Y value (Fmin)- By extrapolating results, we can conclude 
that in these cases helium abundances are lower than y,nin, and 
then obtain the other FSPs. In those cases, however, these stars 
must be brighter or cooler to avoid being on the wrong side of 
the ZAMS for the corresponding Y value. 

Since stars with extrapolated results are found on the fainter 
and hotter side of the ZAMS for the given chemical composition. 
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Fig. 2. Helium abundance Y versus metallicity Z (left panels), and 
stellar mass versus Age (right panels) for stars listed in C06. Circles 
and triangles represent stars with ages younger and older than 13.5 
Gyr, respectively. The Sun's properties (Z = 0.0167, Y = 0.262, 
Age = 4.6 Gyr, and M = I Mq) are represented by the symbol of the 
Sun. Open and filled symbols are the extrapolated and interpolated re- 
sults, respectively. The dotted line represents the reference Y-Z relation 
Y = 0.240 -I- 2.0 X Z. Upper panels: FSPs determined using the method 
from C07. Middle panels: FSPs determined using the method pro- 
posed in this paper. Lower panels: Similar to the middle panels, except 
theoretical evolutionary tracks have a variable [o'/Fe] value (see text for 
more details). 



the main findings exclude all stars with at least one FSP outside 
the respective range. This means that only stars with 0.23 < Y < 
0.37, 0.5 < M/Mq < 1.1, and 1.6 x 10"^ < Z < 6.0 x lO^^ are 
accepted. In other words, stars with open symbols in the middle 
and bottom panels of Fig. |2]are not included in the results of the 
next section (Figs. [3]and|4]i. 



4. Application 

To apply the method proposed in this paper, one additional 
database is selected with different properties to those of C06. 
The second database was created bv lBaumann et al.l (1201 Oi here- 
after BIO), and consists of 117 stars with near-solar proper- 
ties: -0.4 < [Fe/H] < +0.3, 5600 < T^ff[K] < 6100, and 
4.0 < log g < 4.6. These stars have excellent spectroscopic mea- 
surements, where mean errors for the spectroscopic variables are 
smaller than in the C06 catalog, namely cr[pe/H] ~ ±0.025 dex, 
CTj-.ij X +40 K, and o-\osg ~ +0.06 dex. The database contains 
a good set of well-measured stars for comparison with the Sun. 
For each star, the Hipparcos parallax and V magnitude are ob- 
tained from the SIMBAD database, in order to use My as op- 
posed to Mboi. In this case, theoretical models are transformed 
to the observational plane using the bolometric corrections from 
ICasteni&Kuruczl(l2003h . 
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Fig. 3. FSPs of stars from C06 and BIO. The error bars for K, M, and Age are determined based primarily on the error in the surface gravity 
(c\ogg)- The error bars in Z are established by propagating the error in [Fe/H] listed by each respective author. Blue circles and red triangles 
indicate stars (formally) younger and older than the accepted age of the Universe (~13.5 Gyr), respectively. Dotted lines represent the Y-Z relation 
Y = 0.240 + 2.0 X Z. 



4. ) . Fundamental stellar parameters and their errors 

The FSPs for the two selected databases and their respective er- 
rors are shown in Fig. |3] Before discussing these results, it is 
important to note that only non-extrapolated results are being 
used. As such, FSPs are determined for only 33% and 20% of 
the total number of stars in the C06 and BIO databases, respec- 
tively. Using extrapolated results produces higher percent-ages, 
as occurs in the middle panels of Fig. |2] where FSPs are deter- 
mined for 75% of C06 stars. Similarly, the percentage of stars 
younger than 13.5 Gyr with determined results is 25% for C06, 
and 87% for BIO. 

Formal error bars in Y, M, and Age are established based 
on the error in the spectroscopic value of log g combined with 
the same analysis as in the bottom-right panel of Fig. [T] In 
the case of Z, error bars are determined by propagating directly 
from the error in [Fe/H]. Although we have argued against using 
extrapolated results, a linear extrapolation algorithm was applied 
to estimate error bars outside the limits of the theoretical models 
to make sure that errors do not seem smaller than they really are. 

For the selected databases, Y-Z, M-Age, and Z-Age rela- 
tionships are shown in the left, middle and right-hand panels of 
Fig. [3] respectively. Based on the C06 results (upper panels) one 
finds that the reference Y-Z relation (dotted line) seems to fit the 
lower values of Y for every Z value, whereas the maximum Y 



value shows a steeper increase with Z. In addition, no relation- 
ship is observed between Z and Age. However, as previously 
mentioned, errors may occur in calculating Y and M values for 
stars (formally) older than the age of the Universe (red triangles) 
because of difficulties modeling the stellar radius. In principle, 
one can consider such a (formally) very old star as having un- 
derestimated g. Consequently, its helium abundance would be 
overestimated and mass and age underestimated. On the other 
hand, underestimation of its brightness (owing to parallax or the 
brightness itself) or overestimation of T^ff may cause overesti- 
mation of Y and mass, as well as underestimation of age. As 
such, depending on the stellar evolutionary phase of the star and 
its mass the errors in the FSPs may be large or small. 



For the BIO database (lower panels in Fig. |3), a substantial 
percentage of stars have similar metallicities, ages, and masses. 
However, the spread in helium abundance is greater than ex- 
pected when compared to the reference Y-Z relation. Almost 
all stars in this database with established FSPs are younger than 
13.5 Gyr, except for three stars (or one if the error bars are con- 
sidered). This illustrates the good agreement between theory and 
recent observations for stars with masses greater than 0.60 Mq. 
For this database no relationship is found between Z and Age. 
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In Fig. |4] the masses and ages of stars determined in this 
study are compared to those determined by BlCfl For BIO a 7- 
Z rela tion with AF/AZ = 2.0 is assumed for all stars (lYi et al.l 
I2003h . This comparison demonstrates the consequences of de- 
termining masses and ages when assuming that all stars have the 
same Y value for a given Z value. Thus, we can conclude that if 
a universal Y-Z relation is assumed valid for all stars, it is possi- 
ble to identify differences of |AM| a; 0.2 Mq and |AAge| ^ 2 Gyr 
with respect to the case when Y is variable. 

4.2. Is AY/AZ constant? 

As pointed out in Sect. [T] the Y-Z relation is mostly reg- 
ulated by the constant value of AY/AZ. However, several 
authors have emphasized the possibility that different stel- 
lar populations might have diffe rent Y-Z relations, implying 
that AF/AZ is not constant (e.g., Catelan & de Freitas Pachecd 



T996tlCatelanl2009l:lNataf & Udalskil201 UlNataf et al.l2011allbt 



Nataf & Gouldl l2012h . To solve this conundrum, an attempt 
should be made to determine the Y value for each star individ- 
ually, as performed in this study. However, as shown in Fig. |3j 
this method requires highly accurate measurements of surface 
gravity, with o-\o„g < 0.05 dex, otherwise large error bars may 
make it impossible to obtain a reliable conclusion. In the near 
future, high-precision observations and more accurate theoreti- 
cal models could establish whether AY/AZ is not constant only 
in some GCs, or if this is also valid among the field populations 
in galaxies. 

4.3. Consequences of non-canonical effects 

Given that the surface gravity is used here as one of the main 
parameters to determine the FSPs, it is important to know which 
stellar pr operties may pote ntially affect the measurement of g. 
Recentlv lBasu et aTl (1201 2l hereafter B 12) have tested how prop- 
erties of single stars (M, R, and log g) are affected by the uncer- 
tainties in stellar models. They studied the effects induced by 
different metallicity scales, atmospheric models {T -t relation), 
mixing length parameters o-/, and shifts in the temperature scale 
for observed stars with and without seismic data. We focus on 
the last, since our stars also lack seismic data. In this case, their 
results show that the errors in the determination of M and R are 
around 8% and 14%, respectively, for errors in the observational 
data of cTj-^g- - 50 K, o-[fs./u] =0.1 dex, and (Xiog j = 0.1 dex. 

While adopting variable a/ values, as suggested by B12, rep- 
resents an interesting means of improving the estimations of the 
stellar parameters (although in this case, care should be taken to 
properly calibrate the temperature scale implied by the models, 
so that it retains some predictive power for stars other than the 
Sun, which at pres ent remains th e sole calibrator that is used to 
obtain ai; see, e. g.. ICatelanll20 1 2i and references therein), a vari- 
able initial helium abundance represents another potentially fea- 
sible alternative, especially at the main sequence and sub giant 
branch phases. Regarding this topic and because here and in B 12 
non-rotating stellar models are used, one can specul ate that some 
calibration parameters of the model (such as ai, iKapvla et al.l 



* Although, in principle, it is not possible to directly compare the 
results of this study (M and Age) with the respective values deter- 
mined by B IO since their masses and ages are obtained using the 
isochirones jYi et al.l l2001'). substantial differences are not expected, 
in view of the reasonable level o f agreement between Y" and BaSTI 
found inlPietrinfemi et al] ( I2004I) . and between BaSTI and PGPUC in 
IValcarceetal.l(l2012h . 
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Fig. 4. Comparison between masses (upper panel) and ages (lower 
panel), determined using theoretical models that follow a Y-Z relation 
(BIO) and those with non-fixed helium abundance (this work, TW). 
Dotted lines indicate a ratio of 1:1. Below each figure is a histogram 
depicting the absolute differences between both approaches. 



I2OOI may indeed not be valid for all stars, since stellar rotation 
may differ from the rotation pattern which is observed in the Sun 
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(e.g.. lDe Medeiros et al.ll2006t ICortes et al.ll2009h - although in 
practice the deviations from the canonical case may be small. 

The intere s ted re ader is also referred to Sect. 5 in 
iGennaro et al ] (120101) for a detailed discussion of possible 
sources of uncertainties that may be present in comparisons of 
this kind. 

Since one of the aims of this paper is to determine whether 
the Y-Z relation can be used with stellar models with the same 
characteristics as those created by the PGPUC SEC, an analysis 
of these possible differences is beyond the scope of this paper 
Nevertheless, we emphasize the importance of taking Y into ac- 
count as an unknown parameter, when attempting to derive the 
FSPs of stars. 

5. Conclusions 

In this paper a new method for determing the FSPs (7, Z, M, 
and Age) of nearby stars is presented, where the required input 
parameters include the chemical composition ([Fe/H], as well as 
[o'/Fe]), Mboi or the absolute magnitude in a given filter (e.g.. 
My), Teif, and g. 

This method uses almost 600 new evolutionary tracks for 7 
masses (0.5 < M/Mq < 1.1), 7 helium abundances (0.230 < F < 
0.370), and 12 metallicities (1.6 x lO""* < Z < 6.0 x 10"^) foj. 
[a/Fe] - 0.0. These new evolutionary tracks are now included 
in the PGPUC online database, making it possible to interpolate 
evolutionary tracks and isochrones with a variable a-element 
distribution (0.0 < [a/Fe] < +0.3). 

This method is tested using the C06 database, concluding 
that an age of 5 Gyr for all nearby stars is not a good approxi- 
mation, since this assumption leads to an increasing underesti- 
mation of Y as metallicity decreases. Moreover, the method is 
only reliable for masses higher than about 0.60 M©, in view of 
problems in determining the radius of stars below that limit. 

The FSPs are determined for 20% of the stars of the BIO 
spectroscopic database. When comparing the masses and ages 
obtained in this study (where Y is not assumed a priori to de- 
pend on Z) with those obtained by BIO (where Y follows a Y-Z 
relation with SYj t^Z - 2.0), average differences in the masses 
and ages of |AM| k. 0.02 Mq and |AAge| k. 2 Gyr are recorded, 
respectively, for the range of masses used here (0.5 < MIMq < 
1.1). 
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